NYU-TH/00/05/05 
TPI-MINN-16/00 
UMN-TH-1850/00 
hep-th/ 00073451 



Vacuum Structure and the Axion Walls in Gluodynamics 
and QCD with Light Quarks 



Gregory Gabadadze^* and M. Shifman^ 

^Department of Physics, New York University, New York, NY 10003 
"^Theoretical Physics Institute, University of Minnesota, Minneapolis, MN 55455 



Abstract 

Large pure gluodynamics was shown to have a set of metastable vacua with 
the gluonic domain walls interpolating between them. The walls may separate the 
genuine vacuum from an excited one, or two excited vacua which are unstable at 
finite Nc- One may attempt to stabilize them by switching on the axion field. We 
study how the light quarks and the axion affect the structure of the domain walls. 
In pure gluodynamics (with the axion field) the axion walls acquire a very hard 
gluonic core. Thus, we deal with a wall "sandwich" which is stable at finite Nc- In 
the case of the minimal axion, the wall "sandwich" is in fact a "27r" wall, i.e., the 
corresponding field configuration interpolates between identical hadronic vacua. The 
same properties hold in QCD with three light quarks and very large Nc- However, 
in the realistic case of three-color QCD the phase corresponding to the axion field 
profile in the axion wall is screened by a dynamical phase associated with the r]', so 
that the gluon component of the wall is not excited. We propose a toy Lagrangian 
which models these properties and allows one to get exact solutions for the domain 
walls. 
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1 Introduction 



The early studies of the chiral Ward identities in QCD revealed that the vacuum 
energy density depends on the vacuum angle 9 through the ratio 6/Nf, where Nf is 
the number of quarks with mass nig <^ A. Shortly after, Witten and Di Vecchia 
and Veneziano showed that this structure occurs naturally, provided that there 
exist Nf states in the theory such that one of them is the true vacuum, while others 
are local extrema; all are intertwined in the process of "the 6 evolution." Namely, 
in passage from 6 = to 9 = 2it, from 6 = 2-71 to 6 = Att, and so on, the roles of 
the above states interchange: one of the local extrema becomes the global minimum 
and vice versa. This would imply, with necessity, that at 9 = kn (where k is an odd 
integer) there are two degenerate vacuum states. Such a group of intertwined states 
will be referred to as the "vacuum family." The crossover at 9 = n, 3tt, etc. is called 
the Dashen phenomenon 0]. 

This picture was confirmed by a detailed examination of effective chiral La- 
grangians [0, ^, |], ^ (for a recent update see For two and three light quarks 
with equal masses it was found that the vacuum family consists of two or three 
states respectively; one of them is a global minimum of the potential, while others 
are local extrema.Q At 9 = n the levels intersect. Thus, Crewther's dependence 
on 9/Nf emerges. 

On the other hand, the examination of the effective chiral Lagrangian with the 
realistic values of the quark masses, md/rriu ~ 1.8 , 171^/771^ ~ 20, yields 0, ^ [7[ a 
drastically different picture - the vacuum family disappears (shrinks to one state); 
the crossover phenomenon at = vr is gone as well. 

This issue remained in a dormant state for some time. Recently arguments were 
given that the "quasivacua" (i.e. local minima of the energy functional), which to- 
gether with the true vacuum form a vacuum family, is an indispensable feature of 
gluodynamics. The first argument in favor of this picture derives from super- 
symmetric gluodynamics, with supersymmetry softly broken by a gluino mass term. 
The same conclusion was reached in Ref. based on a D-brane construction in the 
limit of large Nc- In fact, one can see that in both approaches the number of states 
in the vacuum family scales as Nc- Finally, an additional argument may be found 
in a cusp structure which develops once one sums up [jlOl sub leading in l/N^ terms 
in the effective rj' Lagrangian. At Nc = 00 the states from the vacuum family are 
stable, and so are the domain walls interpolating between them |T^. 

When Nc < 00 the degeneracy and the vacuum stability is gone, strictly speak- 
ing. It is natural to ask what happens if one switches on the axion field. This 
generically leads to the formation of the axion domain walls. The axion domain 
wall |12[ presents an excellent set-up for studying the properties of the QCD vac- 
uum under the 9 evolution. Indeed, inside the axion wall, the axion field (which, in 

^We stress that the states from the vacuum family need not necessarily lie at the minima of the 
energy functional. As was shown by Smilga |Q, at certain values of 9 some may be maxima. Those 
which intersect at 6 — kn {k odd) are certainly the minima at least in the vicinity of 9 = kir. 
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fact, coincides with an effective 6) changes slowly from zero to 2tt. The character- 
istic length scale, determined by the inverse axion mass m~^, is huge in the units 
of QCD, A~^. Therefore, by visualizing a set of spatial slices parallel to the axion 
wall, separated by distances ^ A^^, one obtains a chain of QCD laboratories with 
distinct values of 6es slowly varying from one slice to another. In the middle of the 

wall 9eS = TT. 

Intuitively, it seems clear that in the middle of the axion wall, the effective value 
of ^eff = TT. Thus, in the central part of the wall the hadronic sector is effectively in 
the regime with two degenerate vacua, which entails a stable gluonic wall as a core 
of the axion wall. In fact, we deal here with an axion wall "sandwich." Its core is 



the so-called D wall, see 13 



Below we will investigate this idea more thoroughly. We also address the question 
whether this phenomenon persists in the theory with light quarks, i.e., in real QCD. 
Certainly, in the limit Nc = oo the presence of quarks is unimportant, and the 
axion wall will continue to contain the D-wall core. As we lower the number of 
colors, however, below some critical number it is inevitable that the regime must 
change, the gluonic core must disappear as a result of the absence of the crossover. 
The parameter governing the change of the regimes is A/N^ as compared to the 
quark mass m^. At rriq <^ A/Nc, even if one forces the axion field to form a wall, 
effectively it is screened by a dynamical phase whose origin can be traced to the 77', 
so that in the central part of the axion wall the hadronic sector does not develop 
two degenerate vacua. The D walls cannot be accessed in this case via the axion 
wall. 

A part of this paper is of a review character. We collect relevant assertions 
scattered in the literature. The main original results - the occurrence of the D-wall 
core inside the axion wall in pure gluodynamics and in QCD with A ^ rrig ^ A/A^c 
- are presented in Sees. 4-7. 

Recently, the issue of hadronic components of the axion wall in the context of 



a potential with cusps |10| was addressed in |14, O, 16|. However, the gluonic 



component of the axion walls was not studied. The t]' component in the axion walls 



was considered in |17, 14]. As far as we understand, in actuality the rj' component 



is totally unstable, and cannot be discussed in the static regime. 



2 Invisible Axion and Axion Walls 

In this section we briefly review the axion set-up, mainly with the purpose of setting 
the relevant notation. 

The axion was originally introduced by Weinberg []TB|| and Wilczek [|TP[ to solve 



the strong CP problem which arises in QCD if physics at very short distances (say, 
of order the Planck scale) generates a non-vanishing 6 term. In the original version 
the axion was coupled directly to the light u, d and s quarks. 

Shortly after, it became clear that the original construction of Weinberg and 
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Wilczek is not viable from the phenomenological standpoint and the axion mech- 
anism was further developed: "invisible axions" were introduced. The concrete 
version we will keep in mind is the KSVZ axion I^U], |2T[| (although other versions 
can be considered too [^]). One introduces a complex scalar field coupled to a 
hypothetical quark field Q in the fundamental representation of the color SU(3), 
with no weak interactions, 

A/: = (t>QRQL + H.c. . (1) 

The modulus of is assumed to develop a large vacuum expectation value //-\/2, 
while the argument of becomes the axion field a, modulo normalization, 

a(x) = fa{x) , a{x) = Arg0(a;) , / » A . (2) 

Then the low-energy coupling of the axion to the gluon field is 



1 ^ 



so that the QCD Lagrangian depends on the combination 6 + a{x). 

In general, one could introduce more than one fundamental field Q, or introduce 
them in a higher representation of the color group. Then, the axion-gluon coupling 
(H) acquires an integer multiplier iV, 



/^C' = -aN^^Gl^Gl^. (4) 



1 N 3' 
— ai\ „ 

/ 327r2 

This is sometimes referred to as the axion index, not to be confused with A/" of 
extended super symmetry, nor with N^, the number of colors. The minimal axion 
corresponds to A^ = 1. In the general case the QCD Lagrangian depends on the 
combination 6 + Na{x). The phenomenon of formation of the axion domain walls 



is being discussed in the literature for a long time |jT2[. The character of the axion 
walls depends on A^. For A^ = 1 there is no physical vacuum degeneracy (except 
aX 6 = tt). Since the wall interpolates between the vacuum and its "27r copy" it 



can be bounded by a closed axion string (see Ref. PB[ for a review). Thus, such 



a wall can have a finite longitudinal extent. This wall is classically unstable as it 
shrinks its size down by emitting axions. Moreover, the A^ = 1 axion walls are, 
strictly speaking, unstable even if they have an infinite extent. They can decay 
quantum mechanically. The decay process is due to tunneling between the identical 
vacua separated by a barrier. In fact, a hole can be created in the wall - a domain 
where the modulus of the field = f /V^ vanishes, and its phase can be "unwind." 
This hole then expands to infinity removing the wall completely. Numerically this 
process is extremely suppressed due to the fact that / is very large in the vacuum, 
and suppressing |0| to make a hole in the wall costs a lot of energy. The suppression 



factor for tunneling was estimated to be ~ exp {(—const / m„ ) ln(/ m^ )} 



Thus, the infinite-extent wall can be considered stable for all practical purposes. 
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If > 2, there is a residual vacuum degeneracy of the type; the walls 
connecting distinct vacua must have infinite area and must be perfectly topologi- 
cally stable (they are cosmologically unacceptable, since they would over-close the 
Universe p^). 

Since we have little to add on the process of the wall formation in the early 
universe, for our purposes - consideration of the walls in the static environment - 
the distinction between = 1 and > 2 is unimportant. For simplicity we will 
deal with the N = 1 axions. All formulae are readily adjustable for N = 2 and 
higher. 



3 Two Scenarios (A Signature of the Hadronic 
Core) 

The invisible axion is very light. Integrating out all other degrees of freedom and 
studying the low-energy axion effective Lagrangian must be a good approximation. 
The axion effective potential in QCD can be of two distinct types. 

Assuming that for all values of 6 the QCD vacuum is unique one arrives at the 
axion effective Lagrangian of the form 

Ca = f [lid.aY + ml (cos(a + 9)- 1)] . (5) 

The axion potential does not have to be (and generically is not) a pure cosine; it 
may have higher harmonics. In the general case it is a smooth periodic function 
of a + ^, with the period 27r. For illustration we presented the potential as a pure 
cosine. This does not change the overall picture in the qualitative aspect. 

As we will see below, a smooth effective potential of the type (|^) emerges even if 
the (hadronic) vacuum family is non-trivial, but the transition between the distinct 
hadronic vacua does not occur inside the axion wall. This is the case with very light 
quarks, -C A/Nc- In the opposite limit, one arrives at the axion potential with 
cusps, considered below. 

In the theory (^ one finds the axion walls interpolating between the vacuum 
state at a = — 6* and the same vacuum state at a = — + 27r, 

a{z) + ^ = 4arctan(e'"'"^) , (6) 

where the wall is assumed to lie in the xy plane, so that the wall profile depends 
only on z. This is the most primitive "27r wall." 

The tension of this wall is obviously of the order of 

Ti~/^m„. (7) 

Taking into account that /^m^ ~ x where x is the topological susceptibility of the 
QCD vacuum, we get 

Ti~xMa- (8) 
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The inverse proportionality to is due to the fact that the transverse size of the 
axion wall is very large. 

Let us now discuss the axion effective potential of the second type. In this case 
the potential has cusps, as is the case in pure gluodynamics, where the axion effective 
Lagrangian is of the form 

C = ^(9,a)^ + min|iVfA^cos^±^} , (9) 

(for a more detailed discussion see below). Here the 9 angle was absorbed in the 
definition of the axion field. The axion wall interpolates between a = and a = 2tt. 

What is the origin of this cusp? The cusps reflect a restructuring in the hadronic 
sector. When one (adiabatically) interpolates in a from to 27r a gluonic order 
parameter, for instance (GG), necessarily experiences a restructuring in the middle 
of the wall corresponding to the restructuring of heavy gluonic degrees of freedom. 
In other words, one jumps from the hadronic vacuum which initially (at a = 0) had 
(GG) = into the vacuum in which initially (GG) ^ 0. Upon arrival to a = 27r, 
we find (GG) = again. This implies that the central part of such an axion wall 
is dominated by a gluonic wall. Thus, the cusp at a = vr generically indicates the 
formation of a hadronic core, the D wall in the case at hand. 

Returning to the question of the tension we note that 

X ~ A^'A^'^ , rria ~ Is^N^f"^ in pure gluodynamics , 
X ~ A^A^c^g, rua ~ A3/2^/2jYi/2y-i in QCD with light quarks, (10) 
which implies, in turn, 

{fA^N^ in pure gluodynamics 
(11) 
jA3/2^i/27vV2 in QCD with hght quarks . 

Here is the light quark mass. 

The presence of the large parameter / in Ti makes the axion halo the domi- 
nant contributor to the wall tension. The contribution of the hadronic component 
contains only hadronic parameters, although it may have a stronger dependence on 
Nc- Examining the cusp with an appropriately high resolution one would observe 
that it is smoothed on the hadronic scale, where the hadronic component of the 
axion wall "sandwich" would become visible. The cusp carries a finite contribution 
to the wall tension which cannot be calculated in the low-energy approximation 



25| . To this end one needs to consider the hadronic core explicitly. The tension 
of the core Tcoro ~ A^A^ci while the tension of the axion halo Thaio ~ /A^ (in pure 
gluodynamics). 

We pause here to make a comment on the literature. The consideration of the 
axion walls in conjunction with hadrons dates back to the work of Huang and Sikivie, 
see Ref. |]12|. This work treats the Weinberg- Wilczek N = 2 axion in QCD with 
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two light flavors, which is replaced by a chiral Lagrangian for the pions, to the 
leading order (quadratic in derivatives and linear in the light quark masses). It is 
well-known [0, |^ that in this theory the crossover phenomenon takes place at 
rriu = rrid- In the realistic situation, (m^ — rriu) / {rrid + rriu) ~ 0.3 considered in Ref. 
|T^, there is no crossover. The pions can be integrated over, leaving one with an 
effective Lagrangian for the axion of the type (|^) (with a — >■ 2a). The potential 
is not pure cosine, higher harmonics occur too. The axion halo exhausts the wall, 
there is no hadronic core in this case. 

At the same time, Huang and Sikivie (see Ref. found an explicit solution for 
the "tt"" component of the wall. In fact, this is an illusion. The Huang-Sikivie (HS) 
solution refers to the bare ir^ field. To find the physical 7r° field one must diagonalize 
the mass matrix at every given value of a (the bare fa is the physical axion field 
up to small corrections ~ f'^l P where /tt stands for the pion decay constant). Once 
this is done, one observes that the physical pion field, which is a combination of the 
bare pion and /a, is not excited in the HS solution. The equation (2.16) in the HS 
paper is exactly the condition of vanishing of the physical pion in the wall profile. 
This explains why the wall thickness in the HS work is of order with no traces 
of the component. The crossover of the hadronic vacua at a = 7r/2 (remember, 
this is = 2 model) could be recovered in the Huang-Sikivie analysis at = m^. 
However, the chiral pion Lagrangian predicts in the two-quark case the vanishing of 
the pion mass in the middle of the wall, for accidental reasons. This is explained in 
detail by A. Smilga, Ref. 0. 



4 Vacuum Structure in Gluodynamics with Invis- 
ible Axion 

First we will summarize arguments in favor of the existence of a nontrivial vacuum 
family in pure gluodynamics. 

The first indication that the crossover phenomenon may exist in gluodynamics 
comes 1^ from supersymmetric Yang-Mills theory, with supersymmetry being broken 
by a gluino mass term. The same conclusion was reached in Ref. based on a 
D-brane construction in the limit of large A^^. In both approaches the number of 
states in the vacuum family is Nc- 

The Lagrangian of softly broken supersymmetric gluodynamics is 

+ ^^^G'^^uG'^^u, (12) 

where m is the gluino mass which is assumed to be small, m <^ A. 

There are distinct chirally asymmetric vacua, which (in the m = limit) are 
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labeled by 



(A2), = iV,A3expfz^^l^V £ = 0,l,...,iV,-l. (13) 



N, 



c 



At m = there are stable domain walls interpolating between them |2^. Setting 



m 7^ we eliminate the vacuum degeneracy. To first order in m the vacuum energy 
density in this theory is 

S = ^(A^) + H.c. = -mAT^A^ cos . (14) 

Degeneracy of the vacua is gone. As a result, all the metastable vacua will decay 
very quickly. Domain walls between them, will be moving toward infinity because 
of the finite energy gradient between two adjacent vacua. Eventually one ends up 
with a single true vacuum state in the whole space. 

For each given value of 6 the ground state energy is given by 

S{e) = min, i^-mN^A^ cos ^^^| . (15) 

At = TT, Svr, we observe the vacuum degeneracy and the crossover phenomenon. 
If there is no phase transition in m, this structure will survive, qualitatively, even 
at large m when the gluinos disappear from the spectrum, and we recover pure 
gluodynamics. 

Based on a D-brane construction Witten showed that in pure SU(A"c) (non- 
supersymmetric) gluodynamics in the limit Nc —>■ oo a. vacuum family does exist :0 
the theory has an infinite group of states (one is the true vacuum, others are non- 
degenerate metastable "vacua") which are intertwined as 6 changes by 27r x (integer), 
with a crossover at ^ = 7rx(odd integer). The energy density of the k-th state from 
the family is 

=iV.'A'f(^), (16) 

where F is some 27r-periodic function, and the truly stable vacuum for each 6 is 
obtained by minimizing Sk with respect to k, 



£{9) = Nl A' miufc F , (17) 




^This was shown in Ref. assuming that there is no phase transition in a certain parameter 
of the corresponding D-brane construction. In terms of gauge theory, this assumption amounts of 
saying that there is no phase transition as one interpolates to the strong couphng constant regime. 
Thus, the arguments of |^ have the same disadvantage as those of SUSY gluodynamics where one 
had to assume the absence of the phase transition in the gluino mass. 
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much in the same way as in Eq. (|I5|). 

At very large Eq. ([T7|) takes the form 



S{e) = minfc {e + + C ( ^ ) • (18) 



The energy density S{6) has its absolute minimum at ^ = 0. At Nc = oo the 
"vacua" belonging to the vacuum family are stable but non-degenerate. To see that 
the lifetime of the metastable "vacuum" goes to infinity in the large limit one 
can consider the domain walls which separate these vacua JIT], These walls are 
seen as wrapped D branes in the construction of [0, and they indeed resemble many 
properties of the QCD D branes on which a QCD string could end. We refer to them 
as D walls because of their striking similarity to D2 branes. The consideration of 
the D walls has been carried out |jll| and leads to the conclusion that the lifetimes 
of the quasivacua go to infinity as exp ( const A^^). 

Moreover, it was argued |13| that the width of these wall scales as 1/Nc both, 
in SUSY and pure gluodynamics. To reconcile this observation with the fact that 
masses of the glueball mesons scale as iV°, we argued |jl3[ that there should exist 



heavy (glue) states with masses oc A^^^ out of which the walls are built. The D-brane 
analysis effective Lagrangian arguments and analysis of the wall junctions pO 



support this interpretation. These heavy states resemble properties of the DO branes. 
The analogy is striking, as the DO branes make the D2 branes from the standpoint 
of the M(atrix) theory [^], so these QCD "zero-branes" make the QCD D2 branes 
(i.e. domain walls). The distinct vacua from the vacuum family differ from each 
other by a restructuring of these heavy degrees of freedom. They are essentially 
decoupled from the glueballs in the large Nc limit. 
Now we switch on the axion 



AC = -f(d,a)(d''a) + — G%,G;, , (19) 



with the purpose of studying the axion walls. The potential energy S{9) in Eq. (|1] 
or (|18D is replaced by £{9 + a). 

Since the hadronic sector exhibits a nontrivial vacuum family and the crossoveiQ 
at 6^ = TT, Stt, etc., strictly speaking, it is impossible to integrate out completely the 
hadronic degrees of freedom in studying the axion walls. If we want to resolve the 
cusp, near the cusp we have to deal with the axion field plus those hadronic degrees 
of freedom which restructure. In the middle of the wall, at a = vr, it is mandatory 
to jump from one hadronic vacuum to another - only then the energy of the overall 
field configuration will be minimized and the wall be stable. Thus, in gluodynamics 
the axion wall acquires a D-wall core by necessity. 

One can still integrate out the heavy degrees of freedom everywhere except a 
narrow strip (of a hadronic size) near the middle of the wall. Assume for simplicity 

■^See also closely related discussions in Ref. 

^For nonminimal axions, with N > 2, the crossover occurs at a = kn/N. 
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that there are two states in the hadronic family. Then the low-energy effective 
Lagrangian for the axion field takes the form (^. The domain wall profile will also 
exhibit a cusp in the second derivative. The wall solution takes the form 

{Sarctan (e™"''tan^) , at z < 
-27r + 8 arctan (e™"^ tan^J , at z > , 

where the wall center is at z = 0. 

Examining this cusp with an appropriately high resolution one would observe 
that it is smoothed on the hadronic scale, where the hadronic component of the 
axion wall "sandwich" would become visible. The cusp carries a finite contribution 
to the wall tension which cannot be calculated in the low-energy approximation but 
can be readily estimated, Tcore ~ A^iVc. 

Below we will examine this core manifestly in a toy solvable model. Before doing 
so, however, we want to elucidate the issue of the peculiar A^^ dependence (or, better 
to say, its absence), in Eq. (|l^). 



5 Description in Terms of a Three-Index Field 

The expression for the vacuum energy density (|TBp seems somewhat puzzling from 
the point of view of the gluon Lagrangian. Indeed, there are — 1 degrees of 
freedom in gluodynamics. Therefore, naively, one expects that the vacuum energy 
density in the large Nc limit scales as ~ N"^. However, the leading term in Eq. (p!8D 
scales as A''^. As a possible explanation, one could think of a colorless massless 
excitation which would give rise to the energy density (0). However, there are no 
physical massless states in gluodynamics. 

The explanation to this apparent puzzle might come if one introduces a colorless 
composite three- index field which does not propagate any physical degrees of freedom 
p3| , On the other hand, this field gives rise |]2^ to precisely the vacuum energy 
(p!8|). In a sense, this field is similar to the photon in (l-l-l)-dimensional QED, where 
a vector particle has no physical degrees of freedom, but it can create a constant 
electric field background which produces a nonzero energy. 

The three-index field in gluodynamics is defined as follows: 

32vr2 '-^Mi^'^M'^ 4! 41 ' 

where H^uap is the field strength for the potential C^ya, and the square brackets 
denote antisymmetrisation over all indices. Hence, the C^ya held can be expressed 
through the gluon fields as follows: 
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Here fabc denote the structure constants of the corresponding gauge group. The 
derivative in this expression is defined as AdB = A{dB) — (dA)B. Note that the 
Cua/B field is not a gauge invariant quantity. If the gauge transformation parameter 
is denoted as A'*, the three-index field transforms as 

CuaP CyaP + duAa/B — daAyp — dfjAau , (23) 

where A^^ oc A'^df3A°' — A'^daA"'. However, the expression for the field strength 
Hfiuais is gauge invariant. 

At energies below A, all massive glueballs decouple from the effective Lagrangian 
of gluodynamics. Thus, no physical excitations are left. However, there should exist 



a kinetic term for the C field in the low-energy Lagrangian This is related 

to the fact that the correlator of the vacuum topological susceptibility x at zero 
momentum is non-zero in gluodynamics. Neglecting all higher derivative terms and 
also terms suppressed in the large Nc limit one arrives at the effective Lagrangian 
for the C field of the form 



_i_ Hi^^, ^ S ^-^^^ . O ij) . (24) 



The first term in this expression reproduces the proper correlation function for the 
topological susceptibility. The second contribution is just the 9 term. Once this 
Lagrangian is set, it is easy to show that the classical equations of motion have a 
constant solution 

H^,uap = - X + 27r/i;) e^^a/B , (25) 

which reproduces the correct large Nc expression for the energy density (0). Note 
that this solution persists even if higher derivatives are included in (^4]). Moreover, 
since the H field does not propagate the dynamical degrees of freedom, the large Nc 
classical solution ( pS]) is also exact quantum-mechanically. 

In this approach, the multiple structure in (^) is related to the quantization 
of the topological charge ||2^. This provides an explanation for the expression (p!8D 
from the point of view of gluodynamics. 

The three-index field C can naturally couple to a D wall. The corresponding 
charge of the D wall is related to the instanton number in gluodynamics ^7j. Thus, 



the D walls are the sources of a constant "electric" field (^) which produces the 
vacuum energy density (0). 

Let us now discuss the mixing of the three-index field with the axion, after the 
latter is switched on. At low energies, when all glueballs are decoupled, two new 
terms emerge in the effective Lagrangian, 



{d,af + ^ ^ . (26) 
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It is known that the pseudoscalar field in four-dimensions is dual to a two-index 
antisymmetric gauge field, B^^ That is to say, the axion Lagrangian (p6D 

can be rewritten in terms of a two-index field. The topological charge density, to 
which the axion is coupled in (^), is rewritten in terms of a three- index field C^j^a 
(prp. It is intriguing to understand what happens with these three- and two-index 
fields after they are coupled to each other (see also a related discussions in \^% ). 
We can rewrite (j2^) in the following equivalent form: 

- \pI + P.d^a + ^ 1^ . (27) 

Here we have introduced an auxiliary field p^. Equations ( p6D and ( pTj) are equivalent 
- to see this one integrates out and substitutes the result p^ = d^a into (^). 

On the other hand, we could first integrate over the axion field in (P?!). This 
gives rise to the following relation: 

1 c"^"^ + 

= y , (28) 

where we have introduced an antisymmetric two- index field Bap- Using the relation 
(p8|), we find that (B^) (or equivalently (p6D) is proportional to 



1 / N 2 

{Cual3 + d[„Ba/3]j ■ (29) 

The sum in the parenthesis is gauge invariant. In fact, it is invariant under both, the 
Abelian transformations on the B field, and the non-Abelian transformations of glu- 
ons. As we mentioned earlier, this latter transformation does not leave C invariant. 
The invariance in (|29| ) is restored, however, due to compensating transformations of 



the B field 38 



At low energies, when all glueballs are decoupled, the expression ( p9D should be 
combined with the 6 term and the gauge invariant kinetic term for C given in (0). 
As a result, the expression (|29| ) is nothing but the gauge invariant mass term for 
the three-index field which is a superposition of C and B fields. In other wards, a 
mixed state of the C field and the B field produces a state with the mass 

= f2- (30) 

This is the physical axion (similar results were first obtained in a different context in 
Ref . [0 by studying correlation functions of the three- index field. This is equivalent 
to the effective Lagrangian approach adopted here). 

Summarizing, we started from gluodynamics where the C field had no physical 
components. The D walls were the sources of the C field. After the axion (rep- 
resented by B) is switched on, the C field and the bare axion mix. The mixed 
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three-index field becomes massive and propagates one massive physical degree of 
freedom, the physical axion. 

The direct physical consequence of this phenomenon is that the three-form charge 
of a D wall in a theory with the axion is screened. As a result, there will be a 
stationary and stable wall in the theory - a superposition of the axion and D wall 
in its core. In the next section we will explicitly find this domain wall "sandwich" 
in a toy model. 



6 An Illustrative Model 

To quantitatively describe the axion walls with the D wall core one has to solve QCD, 
which is way beyond our possibilities. Our task is more modest. We would like to 
obtain a qualitative description of the axion wall sandwich which, with luck, can 
become semi-quantitative. To this end we want to develop toy models. An obvious 
requirement to any toy model is that it must qualitatively reproduce the basic 
features of the vacuum structure which we expect in QCD. In SUSY gluodynamics it 
was possible to write down a toy model with a Z symmetry which "integrates 
in" the heavy degrees of freedom and allows one to investigate the BPS domain 
walls in the large Nc limit p8| (see also |^). We will suggest a similar model in 
(non-supersymmetric) QCD, then switch on axions, and study the axion domain 
walls in a semi-realistic setting. In this model we will be able to find exact solutions 
for the D walls and the axion walls. 

Here is our a simple toy model which has a proper vacuum structure. If an 
appropriate (complex) glue order parameter is denoted by $, the modulus and phase 
of this field describe the 0+"'" and 0^+ channels of the theory, respectively. The toy 
model Lagrangian is 



C 



Ar2(9,,$)*(9,,$)-r($,$* 



Nc-it 



1 + — (l-<l>'^^e-''') 



V = Vq + Vi 



(31) 



Here A is a numerical constant of order one, and x is the vacuum topological sus- 
ceptibility in pure gluodynamics (note that x is independent of N^). The scale 
parameter A is set to unity. 

This model has the vacuum family composed of N^. states. Indeed, the minima 
of the energy are determined from the equations 



dV_ 
9$ 



dV 
(9$* 







(32) 
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which have the following solutions: 



«>£vac = exp (^t j , £ = 0, 1, iV, - 1 . (33) 

In the £-th minimum Vq vanishes, while Vi produces a non- vanishing vacuum energy 
density, 

e,.,N![l-.o.[^)]. (34, 

For each given 6 the genuine vacuum is found by minimization, 

8(6) = AT^xmin, |l - cos } ■ (35) 

The remaining Nc — I minima are quasivacua. Once the heavy field $ is integrated 
out, the vacuum energy is given by the expression (pSf); it has cusps at 6* = vr, 37r 
and so on. Needless to say that the potential (PT| ) has no cusps. 

We will first consider the model (0) without the axion field, at ^ = 0, in the 
limit Nc = oo. In this limit the false vacua from the vacuum family are stable. 

The classical equation of motion defining the wall is 

where primes denote differentiation with respect to z (we look for a solution which 
depends on the z coordinate only). 

This is a differential equation of the second order. It is possible, however, to 
reduce it to a first order equation. Indeed, Eq. ( PE| ) has an obvious "integral of 
motion" ("energy"), 

Ar2 _ y = Const = , (37) 

where the second equality follows from the boundary conditions. In the large 
limit one can parametrize the field $ as follows (p ~ 1): 

$ ^ 1 + ^ . (38) 



Taking the square root of Eq. (|37D , substituting Eq. (|38D and neglecting the terms 
of the sub leading order in 1/Nc we arrive at 

p = iANc (1 - expp) . (39) 



The phase on the right-hand side can be chosen arbitrarily. The choice in Eq. ( pOj) 
is made in such a way as to make it compatible with the boundary conditions for 
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the wall interpolating between $vac = 1 and •I'vac = exp(27ri/Ai'c). This is precisely 
the expression that defines the domain walls in SUSY gluodynamics 0. It is 
not surprising that the same equation determines the D walls in non-SUSY gluody- 
namics - the fermion-induced effects are not important for the D walls in the large 
Nr limit. 



The solution of this equation was obtained in [4^]. In the parametrization p 
a + IT the solution takes the form: 



COST = (cr + 1) exp(— 0" 

ra{z) 
r(0) 



r^'^[exp(2t) - (1 + dt = -ANc\z\ . (40) 

Ja(0] 



The real part of p is a bell-shaped function with an extremum at zero; it vanishes 
at ±oo. The imaginary part of p, on the other hand, changes its value from to 
27T. This determines a D wall in the large Nc gluodynamics. The width of the wall 
scales as l/N^. 

The solution presented above is exactly the same as in SUSY gluodynamics. 
This is not surprising since the ansatz (|38D implies that Vi does not affect the 
solution - its impact is subleading in 1/Nc, while Vq is exactly the same as in the 
SUSY-gluodynamics-inspired model of Ref. ||2^. Moreover, for the same reason the 
domain wall junctions emerging in this model will be exactly the same as in the 



SUSY-gluodynamics-inspired model [jT3[. Inclusion of Vi in the subleading order 
makes the wall to decay. 

Inclusion of the N = 1 axion field amounts to the replacement 

e^e + a 

in Eq. (^Tj), plus the axion kinetic term 

C^n = (^IL±J^^ (d,af + zN,{d,a){^*d,^ - $9^$*) . (41) 

The occurrence of the mixing between a and the phase of $ is necessary, as is readily 
seen from the softly broken SUSY gluodynamics. (To get the potential of the type 
(0) in this model, one must eliminate the GG term by a chiral rotation. Then 
m — i> m exp{{9 + a) / Nc) and, additionally one gets df^ax [the gluino axial current].) 
The term 2$*$ in the brackets has to be included to reproduce the correct mass 
for the axion after the physical heavy state is integrated out. The presence of this 
term signals that QCD dynamics generates not only the potential for the axion but 
also modifies its kinetic term. On the other hand, since $*$ < and, moreover, 
A ^ /, this term can be neglected for all practical purposes. 

We are interested in the configuration with a interpolating between and 27r. 
The phase of $ will first adiabatically follow a/N^, then at a ~ vr, when the phase 
of $ is close to n/Nc, it will very quickly jump by —27c/Nc, and then it will continue 
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to grow as a/Nc, so that when a reaches 27r the phase of $ returns to zero. This 
jump is continuous, although it occurs at a scale much shorter than m~^. This 
imitates the D-wall core of the axion wall. One cannot avoid forming this core, 
since otherwise the interpolation would not connect degenerate states - on one side 
of the wall we would have (hadronic) vacuum, on the other side an excited state. 

In the large Nc limit one can be somewhat more quantitative. Indeed, in this 
approximation the model admits the exact solutions. The gluonic core of the wall 
has the same form as before, Eq. (^), but the phase r is now substituted by the 
superposition t — {a + 6) since the axion field is mixed with the phase of the $ field. 

This very narrow core is surrounded by a diffused axion halo. The axion field is 
described in this halo by the solution to the Lagrangian (^). This takes the form: 

e + a{z) = -27r + 4iVc arctan (e"'"'^ tan-^j , z > 0, 

e + a{z) = - arctan (e'"''^' ^^^^) ' < ^ • (^2) 

Thus, we find explicitly the stable axion wall with a D-wall core. Note that this is 
a usual "27r" wall as it separates two identical hadronic vacua. As we discussed in 
the introduction, this wall is harmless cosmologically. It will be produced bounded 
by global axion strings in the early universe. Bounded walls shrink very quickly by 
decaying into axions and hadrons. 



7 QCD with Three Light Quarks and Axion 

So far we discussed pure gluodynamics with the axion. Our final goal is to study 
QCD with Nf = 3. There are two, physically distinct regimes to be considered in 
this case. In real QCD 

rriu, mrf -C ~ ^ , m„, m^, < A . (43) 

In this regime the consideration of the chiral Lagrangians 0, ^ 0, does not exhibit 
the vacuum family. We will comment on why the light quarks screen the vacuum 
family of the glue sector, so that the axion domain wall provides no access to it. In 
the limit (^) the effects due to the D walls will be marginal. 
On the other hand, in the genuinely large Nc limit 

A 

— < rriu, ma <^ nis <^ A, (44) 

physics is rather similar to that of pure gluodynamics. The light quarks are too 
heavy to screen the vacuum family of the glue sector. 

In what follows we study the axion walls and their hadronic components in the 
limits (|43|) and (H), separately. 
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7.1 One Light Quark 

To warm up, let us start from the theory with one hght quark. In the hmit of large 
Nc this introduces a light meson, "r/"'. An appropriate effective Lagrangian can 
be obtained by combining the vacuum energy density of gluodynamics with what 
remains from the Witten-Di Vecchia-Veneziano Lagrangian at Nj = 1, 



V = -niqA^Nc cos /? + mm < -N^M^ cos — ) . (45) 



' c 



Here /3 is the phase of f/ ~ qi^R-, while ~ K^N^ is the "r^"' coupling constant 
squared. The product Ff3 is the "r^"' field. The first term in V corresponds to the 
quark mass term, JVlU + h.c. (see Eq. (7) in Witten's paper 0). hi Nc = oo 
the second term in V becomes {(3 + 6*)^. It corresponds to (zln dett/ + 9)"^ in Eq. 
(11) in 0. The subleading in 1/A^c terms sum up into a 2ti periodic function of the 
cosine type, with the cusps. It is unimportant that we used cosine in Eq. (|45|). Any 
27r periodic function of this type would lead to the same conclusions. The second 
term in Eq. (^5]) differs from the vacuum energy density in gluodynamics by the 
replacement 9 ^ j3 + 9. 

If rriq ^ ^/Nc, the first term in is a small perturbation; therefore, in the 
vacuum, P + 9 = 2nk, and, hence, the 9 dependence of the vacuum energy is 

£.Ud) = -m,A'NcCos9. (46) 

It is smooth, 2tt periodic and proportional to rrig as it should be on general grounds 
in the theory with one light quark. 

The condition nig <^ A/Nc precludes us from sending Nc oo. The would be 
"27r" wall in the variable f3 is expected to be unstable. This is due to the fact that 
at Nc ^ 3 the absolute value of the quark condensate ipip is not "harder" than the 
phase of the condensate j3, and the barrier preventing the creation of holes in the 
"27r" wall is practically absent. 

If one closes one's eyes on this instability one can estimate that the tension of 
the "77"' wall is proportional to A^Nc^'^, with a small correction nigA'^Nc^'^ from the 
quark mass term. The tension of the D-wall core is, as previously, A^A^^c. 

In the opposite limit 

mg > but nig still < A , (47) 

the situation is trickier. Now the first term in V is dominant, while the second is a 
small perturbation. There are Nc distinct vacua in the theory. 
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Then the 6 dependence of the vacuum energy density is 



S^U9)=A'mm{9 + 2ni)\ (49) 

this is similar to that in the theory without hght quarks (i.e., the same as in glu- 
odynamics). The "r/"' wall is stable at Nc —>■ oo, with a D-wall core in its center. 
The 7]' wall is a "27r" wall. 



From this standpoint, the quark with the mass (^) is already heavy, although 
the is still light on the scale of A, 

M,, Ai/2<< A. 

So far the axion was switched off. What changes if one includes it in the theory? 
The Lagrangian now becomes 

V = -m.A^N, cos (3 + mm | - A^^^A' cos + ^+ | ^ (50) 

where the 6 angle is absorbed in the definition of the axion field. 

The bare "77"' mixes with the bare axion. It is easy to see that in the limit 
rriq <^ A/Nc the physical "77"' is proportional to P + a, rather than to p. Therefore, 
even if we force the axion wall to develop, (i.e. a to evolve from to 2tt) the "77"' 
wall need not develop. It is energetically expedient to have P + a = 0. Thus, the 
effect of the axion field on the hadronic sector is totally screened by a dynamical 
phase P coming from the quark condensate. In other words, the axion wall with the 
lowest tension corresponds to the frozen physical "r^"'. 



P + a = 0. 

There is no hadronic core. The tension of this wall is determined from the term 

OC TTlgA^Nc. 

[If one wishes, one could add an (unstable) "77"' wall to the axion wall. Then 
the "77"' wall, with the D-wall core will appear in the middle of the axion wall, but 
they are basically unrelated. This will be a secondary phenomenon, and the D wall 
core will be, in fact, the core of the "7/"' wall rather than the axion wall.] 

If the quark mass is such that (0) applies, then the axion field a cannot be 
screened, since we cannot freeze P + a everywhere in the axion wall profile at zero - 
at TJiq A/Nc, P is proportional to the physical "77"' and is much heavier than the 
axion field. Thus, in this case the axion wall will be described by the Lagrangian 
(^) and will have a D-wall core. One may also add, on top of it, the "77"' wall. This 
will cost my^A^^'^Nc in the wall tension - still much less than A^A^^c of the D-wall 
core of the axion wall. 
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The limit { ^7\ ) is unrealistic. Moreover, in this limit the D walls taken in isolation, 
without the axion walls, are stable by themselves, although they interpolate between 
nondegenerate states JTT|. 



7.2 Three Light Quarks 

Let us turn to the case of three light flavors. The physical picture is quite similar 
to that of the one-flavor case, see Sec. 7.1. 

We assume the mass matrix Ai in the meson Lagrangian to be diagonal. There- 
fore, we will look for a diagonal U (3) meson matrix which minimizes the potential, 

U = diag (e^'^S e''^\ e'''" ) . (51) 
The potential takes the form 

V = -^m,A=^iVeCos0, + mmj-iVfA^cos 5iA±i±M| . (52) 

As before, we will consider two limiting cases, (|4^) and (^^. 

Let us switch off the axion field first. In the limit of genuinely light quarks, Eq. 
(^), when the second term in the potential (|5^) is dominant, the solutions for (p^s 
were found in [0, They satisfy the relation ^3 ^ and (pi + (p2 = —0- The 
corresponding expression for the vacuum energy density is 

^vac(^) = -N^A^ ^ ml + ml + 2mu m^cos^ . (53) 

As in Sec. 7.1, we deal here with a smooth single- valued function of 6. The inclusion 
of the axion replaces 6 ^ 6 + a a. The physical 77' field is given by the sum 
J2i (pi + ct- It is energetically favorable to freeze this state. Thus, the situation is 
identical to that in the one-flavor case: even if the axion wall is forced to develop, 
the physical t]' wall (which is now the J2i<Pi + Ci wall) does not have to occur. The r]' 
wall is a gateway to the D wall. In the theory at hand the vacuum angle is screened 
in the axion wall, and there is no D-wall core. 

If, nonetheless, the rj' wall is formed due to some cosmological initial conditions, it 
will have a D-wall core (albeit the 1]' wall is unstable in the limit at hand and cannot 
be considered in the static approximation). The would-be t]' wall is independent of 
the axion wall; its effect on the axion wall formation is rather irrelevant. 

In addition to this, a "27?" wall could develop built of nonsinglet mesons, at 
certain values of the quark masses. There is nothing new we could add to this issue 
which is decoupled from the issue of the vacuum family in the glue sector and D 
walls. 

We now pass to the opposite limit (^), when the first term in the potential ( |52[ ) 
is dominant. As in Sec. 7.1, there are distinct vacua with the energy given by 
(^). It is straightforward to show that the potential for the axion in this case is 
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of the form with the cusps which signal the presence of the D-wall core. This 
is similar to what happens in gluodynamics. One cannot avoid having an t]' wall 
in the middle of the axion wall, which entails a D wall too. The D walls separate 
the degenerate vacua. Since they "live" in the middle of the axion wall, they are 
perfectly stable. 

(In addition, there can be "27r" walls in either of 0's or their linear combinations. 
However, these latter are unstable and do not appear in the physical spectrum of 
the theory.) 



8 Conclusions 

Summarizing, we have found that the presence of the axion field and the axion wall 
makes the D wall perfectly stable in gluodynamics at finite Nc. The D wall develops 
as a core of the axion wall. It is unavoidable. 

In QCD with light quarks the axion wall may or may not generate the D-wall 
core. Everything depends on the interplay between the quark masses, A and Nc- In 
the realistic case of genuinely light quarks, see Eq. (^31), the phase associated with 
the axion field in the wall profile is screened by a dynamical phase (which can be 
traced back to the presence of rj'). The rj' is not excited, and neither is the D wall. 
There is no D-wall core in the axion wall. 

If Nc is increased so that Eq. ( ^41) holds the picture changes essentially to that 



one deals with in gluodynamics: the rj' wall is excited, opening the access to the D 
wall. The D-wall core develops in the central part of the axion wall. Unfortunately, 
this limit is unrealistic. 
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